Verfahren und Vorrichtung zur automatischen 
Entzerrung von einkanaligen oder mehrkanaligen 

Bild rn 

5 Die Erfindung bezieht sich auf ein Verfahren und eine Vorrichtung zur automatischen 
Entzerrung von einkanaligen oder mehrkanaligen Bildern, wobei ein Bild auf ein 
Referenzbild ohne genaue Kenntnis der Abbildungsfunktion entzerrt wird. 

Eine ungenaue und/oder unzureichende Kenntnis der Abbildungsfunktion liegt 
10 insbesondere dann vor, wenn Bilder mit bewegten Aufnahmesystemen und/oder einem 
bewegten Bildfeld erfasst werden und die Orientierung des Aufnahmesystems 
gegenuber dem Bildfeld nicht oder nur unter groliem Aufwand eindeutig bestimmbar 
ist. So sind beispielsweise Abbildungsfunktion zur Entzerrung von Satellitenbildern fur 
eine Landvermessung nicht oder nur in Grenzen bekannt. Es ist auch denkbar, dass 
15 eine Verzerrung des Referenzbildes nicht oder nur in Grenzen bekannt ist. 

Die Bilder werden durch ein Aufnahmegerat erfasst und liegen fur die 
Weiterverarbeitung als einkanalige Bilder, beispielsweise Grauwertbilder, oder 
mehrkanalige Bilder, beispielsweise Farbbilder, vor. 

20 

Bilder zur Aufnahme und Abbildung der Erdoberflache werden von bewegten 
Aufnahmesystemen aufgenommen. Die Aufnahmesysteme befinden sich 
beispielsweise auf Satelliten oder Flugzeugen. Bekannte Verfahren zur Aufnahme und 
Abbildung der Erdoberflache sind sogenannte photogrammetrische Verfahren, wobei 

25 die aufgenommenen Bilder fur eine Ubertragung in eine zu erstellende Karte entzerrt 
werden. Die notwendige Transformation eines aufgenommenen Bildes in eine 
gewunschte Darstellung wird als Abbildungsvorgang bezeichnet. Fur die Entzerrung 
wird der Abbildungsvorgang als kollineare Gleichung modelliert. Es existiert bereits 
eine Reihe grundlegender Verfahren zur Georeferenzierung von 

30 Fernerkundungsdaten, die entweder auf dem Prinzip der Positionsberechnung mit Hilfe 
eines Modells des Abbildungsvorgangs (Abbildungsmodells) oder auf dem 
Passpunktprinzip mit Hilfe eines Referenzbildes beruhen. Auch eine Verkopplung 
beider Methoden ist moglich. 

35 Fur die Entzerrung unter Verwendung eines Abbildungsmodells wird der 
Abbildungsvorgang als kollineare Gleichung beschrieben oder mittels einfacherer 
Modelle genahert. Sind alle Parameter des Abbildungsvorgangs bekannt, so sind sehr 
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prazise Bilder erzeugbar. In der Praxis sind jedoch nicht alle Parameter, die bei der 
Bildentstehung eine Rolle spielen, wie z.B. Ort und Lage des Satelliten, ausreichend 
genau bekannt. Ursachlich hierfur ist eine Beschrankung der Messgenauigkeit sowie 
die Tatsache, dass nicht alle Parameter zeitlich stabil sind. Die zeitliche Instability 
5 ergibt sich beispielsweise aus Schwankungen der Satellitenbahn. Weiter erfolgt eine 
Vermessung der Satellitenbahn diskret, wobei die Interpolation zwischen diesen 
Zeitpunkten zusatzliche Fehler ergibt. 

Daneben ist es bekannt, bei unvollstandigen oder ungenauen Abbildungsparametern 

10 eine Entzerrung mit Hilfe von Passpunkten durchzufuhren. Das Verfahren wird 
insbesondere dann eingesetzt, wenn das zu entzerrende Bild, beispielsweise ein 
Satellitenbild, nicht auf eine definierte Geometrie, beispielsweise eine Karte, sondern 
auf ein Referenzbild mit unbekannter Verzerrung abgebildet wird. Passpunkte sind 
kontrastreiche Bildstrukturen, welche zeitlich stabil sind und eine Ortsbestimmung 

15 mittels Korrelation ermoglichen. Die Interpolation zwischen den Passpunkten wird in 
der Regel mit zweidimensionalen Polynomfunktionen vorgenommen. Eine Fehlerquelle 
bei der Passpunktentzerrung ist die Unahnlichkeit der Passpunktstrukturen in den 
verschiedenen Bildern, die hervorgerufen werden durch unterschiedliche 
Aufnahmegeometrien, einen unterschiedlichen Zustand der Objekte, welche die 

20 Passpunktstruktur bilden, oder durch Mittelung unterschiedlicher Objekte oder 
Objektteile. Eine weitere Fehlerquelle ist die fiir die Interpolation eingesetzte 
Polynomfunktion, durch welche in der Regel der Abbildungsvorgang nicht ausreichend 
genau modellierbar ist. Dies gilt insbesondere bei grofien Flachen. Je weiter ein Punkt 
von einem Passpunkt entfernt liegt, desto grofier sind die Auswirkungen aufgrund von 

25 Ungenauigkeiten in der Polynomfunktion. Dieser Effekt zeigt sich vor allem bei 
Polynomen hoheren Grades. 

Aufierdem ist es bekannt, die beiden Methoden, die parametrische Entzerrung und die 
Entzerrung mittels Polynomfunktionen zu verknupfen. 

30 

Ein Nachteil der bekannten Verfahren liegt in einer schlechten Automatisierbarkeit, 
insbesondere bei einem geringen Kenntnisstand des Abbildungsvorgangs. Eine 
Anpassung aller Abbildungsparameter an eine aktuelle Situation ist hierbei meistens zu 
aufwendig, so dass eine reine P olynomentzerrung eingesetzt wird. Die notwendigen 
35 Passpunkte werden dabei einer Datenbank entnommen. Die bisher bekannten 
Verfahren weisen Restfehler in der Positionsbestimmung der Passpunkte des zu 
bearbeitenden Datensatzes auf. Dieser Restfehler lasst sich nicht oder nur sehr 
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aufwendig analytisch beschreiben, so dass auch eine entsprechende Anpassung durch 
Korrektur des Abbildungsmodells nicht moglich ist. 

Der Erfindung liegt daher das technische Problem zugrunde, ein Verfahren und eine 
5 Vorrichtung zur automatischen Entzerrung von Bildern zu schaffen, durch welche ohne 
genaue Kenntnis des Abbildungsmodells ein verbleibender Restfehler verringert wird. 

Die Losung des Problems ergibt sich durch ein Verfahren mit den Merkmalen der 
Anspruche 1 und 9. Weitere vorteilhafte Ausgestaltungen der Erfindung ergeben sich 
10 aus den Unteransprtichen. 

Zur automatischen Entzerrung eines Bild auf ein Referenzbild, wobei Parameter einer 
Abbildungsfunktion der Entzerrung mindestens teilweise unbekannt sind, werden aus 
dem Bild mindestens drei Objekte extrahiert und mindestens drei Passpunkte ermittelt, 

15 wobei charakteristische Objektpunkte der extrahierten Objekte als Passpunkte ermittelt 
werden. Die extrahierten Objekte werden Objekten im Referenzbild zugeordnet, wobei 
eine Zuordnung von Objekten in beiden Bildern aufgrund einer Ahnlichkeit von 
Objekten und/oder eines Vektornetzes erfolgt und das Vektornetz durch eine 
Verbindung charakteristischer Objektpunkte gebildet wird. Anhand der Passpunkte 

20 sowie den entsprechenden Punkten im Referenzbild werden die Parameter der 
Abbildungsfunktion angepasst und/oder es wird eine geeignete bekannte 
Abbildungsfunktion ermittelt. Die Anpassung und/oder Ermittlung erfolgt d erart, dass 
ein Fehlersummenminimum in Positionsdifferenzen zwischen projizierten Passpunkten 
und entsprechenden Passpunkten im Referenzbild erreicht wird. Abbildungsfunktionen 

25 mit festen Parametersatzen sind beispielsweise in einer Datenbank abgelegt. Eine 
derartige Anpassung ist von Interesse, wenn vielfach gleiche Abbildungsvorgange 
vorliegen. Fur eine Wahl geeigneter Objekte ist vorteilhaft, wenn die Objekte markante 
und einmalige Formen aufweisen, so dass grobe Fehler ausgeschlossen werden. Das 
Referenzbild kann dabei entweder in einer bekannten Kartenprojektion abgebildet sein 

30 oder eine unbekannte Abbildungsgrundlage aufweisen. Ist die Abbildungsgrundlage 
bekannt, so ist die Abbildungsfunktion aus einem Abbildungsmodell gegeben. Bei einer 
unbekannten Abbildungsgrundlage des Referenzbildes ist die Abbildungsfunktion eine 
affine Transformation. D urch d ie A npassung d er P arameter u nter B erechnung e ines 
Fehlersummenminimums ist eine Anpassung der Abbildungsfunktion ohne genaue 

35 Kenntnis des Abbildungsmodells und/oder im Falle der fehlenden Abbildungsgrundlage 
ohne Abbildungsmodell moglich. 
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In einer bevorzugten Ausfuhrungsform werden in einem weiteren Verarbeitungsschritt 
die fur eine Anpassung der Abbildungsfunktion verwendeten Passpunkte gewichtet 
Hierfur werden zunachst urn einen moglichen Passpunkt des zu entzerrenden Bildes 
und/oder des Referenzbildes Flachen gebildet, sogenannte Passpunktstrukturen. Die 
5 Flachen sind in der Regel quadratisch mit N*N Pixel, wobei N eine beliebige naturliche 
Zahl ist. Die Passpunktstrukturen werden in das jeweils andere Bild mittels der 
Abbildungsfunktion projiziert, wobei eine Qualitat der Passpunktstruktur mindestens 
durch Veranderlichkeit, Richtungskontrast und/oder Ahnlichkeit beschrieben ist. Die 
Veranderlichkeit wird aus der bereits bestimmten Klassenzugehorigkeit abgeleitet 

10 Klassen mit geringer Veranderlichkeit, beispielsweise Wasser, erzeugen eine hohe 
Qualitat. Objekte mit starker Veranderlichkeit, beispielsweise Vegetationen, haben 
hingegen eine geringe Qualitat. Die Werte fur eine derartige Qualitatsbewertung lassen 
sich vorab bestimmen und in Tabellen ablegen. Ein weiteres Kriterium fur die Qualitat 
eines Passpunktes ist der Richtungskontrast. Fur eine exakte Bestimmung von 

15 Positionen der Passpunkte ist ein grofier Kontrast innerhalb der Passpunktstruktur 
erforderlich. Strukturen mit grolien Kontrasten in horizontaler (X-) Richtung und 
vertikaler (Y-) Richtung der Bildebene haben somit eine hohe Qualitat. Der Kontrast 
wird aus maximalen Grauwertdifferenzen in den Zeilen oder Spalten errechnet. 
Grauwerte werden auf einem oder mehreren Kanalen erfasst, beispielsweise auf drei 

20 Kanalen als so genannte Farbwerte. Bei einer Betrachtung der G rauwertdifferenzen 
auf einem Kanal kann entweder der Kontrast aus dem Bild oder bevorzugt aus dem 
Bild sowie dem zugehorigen Referenzbild in das Qualitatsmafi eingehen. Zudem ist die 
Ahnlichkeit Kriterium fur das QualitatsmafJ. Die Ahnlichkeit wird aus einem Vergleich 
der Form des Grauwertverlaufes in X- und Y-Richtung innerhalb der 

25 Passpunktstrukturen des Grauwertbildes und des Referenzbildes abgeleitet. Wahrend 
Veranderlichkeit und Richtungskontrast im zu entzerrenden Bild und/oder im 
Referenzbild festgelegt sind, lasst sich die Ahnlichkeit durch eine Anpassung 
verbessern. Die einzelnen Qualitatsmafte werden zu einer gemeinsamen 
Qualitatsbewertung verknupft. Fur eine derartige Verknupfung ist eine Vielzahl an, 

30 Methoden denkbar, wobei die einzelnen Qualitatsmafie je nach Aufgabestellung 
unterschiedlich gewichtet werden konnen. Aufgrund der gemeinsamen 
Qualitatsbewertung wird jedem Passpunkt eine Wichtung zugeordnet, mit welcher 
dieser in eine Anpassung der Parameter fur eine Bestimmung der Abbildungsfunktion 
eingeht. Passpunkte, welche ein sehr niedriges Qualitatsmafi aufweisen, sind fur eine 

35 weitere Bearbeitung zu verwerfen. 
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In einer weiteren Ausfuhrungsform erfolgt eine genauere Bestimmung der Position der 
ermittelten Passpunkte im zu entzerrenden Bild durch einen Vergleich der 
Passpunktstrukturen. Dabei ist sowohl eine Projektion von Teilen der 
Passpunktstruktur des Bildes in das Referenzbild, als auch umgekehrt denkbar. 
5 Bevorzugt werden fur beide Bildausschnitte Bildscharfen berechnet und der 
unscharfere Bildausschnitt in den scharferen projiziert. Zunachst erfolgt eine 
Anpassung der Grauwertverteilung oder der Grauwertverteilungen der projizierten 
Bildstruktur an die zugehorige Passpunktstruktur im zu entzerrenden Bild und/oder im 
Referenzbild. Dadurch lassen sich Unterschiede ausgleichen, welche beispielsweise 

10 aufgrund unterschiedlicher Belichtungsverhaltnisse bei der Aufnahme des jeweiligen 
Bildes vorliegen. Anschliefiend werden die Grauwerte der einzelnen Pixel in der 
Passpunktstruktur Oder der Bildstruktur bestimmt und die Differenz zwischen 
benachbarten Pixel gebildet. Diese Differenz wird mit der Differenz der 
korrespondierenden Pixel in der anderen Struktur verglichen und daraus ein Fehler 

15 abgeleitet. Anschliefiend wird die Passpunktstruktur im zu entzerrenden Bild in X- 
und/oder Y-Richtung versetzt und erneut ein Fehler ermittelt. Durch schrittweisen 
Versatz der Passpunktstruktur inX- und Y-Richtung wird eine Fehlermatrix erzeugt. 
Die Schrittgrofie ist dabei geringer als die Ausdehnung der Passpunktstruktur. Als am 
besten angepasste Position wird die Position der Passpunktstruktur identifiziert, welche 

20 den geringsten Fehler aufweist. Positionen, welche auf diese Art ermittelt wurden, 
weisen eine hohe Positionsgenauigkeit auf. Dadurch kann eine Anpassung der 
Abbildungsparameter verbessert werden. 

In einer bevorzugten Ausfuhrungsform wird eine Anpassung der Parameter der 
25 Abbildungsfunktion durchgefuhrt, wobei die Abbildungsfunktion durch eine Anderung 
der Parameter derart veranderbar ist, dass ein Fehlersummenminimum in den 
Positionsdifferenzen zwischen Passpunkten und zugehorigen projizierten Bildpunkten 
unter Berucksichtigung der Wichtung der Passpunkte minimiert wird. Durch eine zwei- 
stufige Anpassung der Abbildungsparameter erfolgt zunachst eine g robe Entzerrung 
30 des Bildes, welche in einem weiteren Schritt verfeinert wird. Es ist jedoch auch 
denkbar, dass Passpunkte bereits in Grenzen bekannt sind und so auf eine erste 
Anpassung verzichtet werden kann. 

In einer weiteren bevorzugten Ausfuhrungsform wird eine Ausgleichsrechnung 
35. durchgefuhrt, wobei der Restfehler in der Positionsbestimmung weiter verringert wird. 
Hierbei wird fur jeden Passpunkt mindestens ein vertikaler und/oder ein horizontaler 
Korrekturwert bestimmt, welcher die Abweichung des Wertes der Abbildungsfunktion 
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vom Wert der Ausgleichsfunktion am Ort des zugehorigen Bildpunktes bei einer 
Projektion mit der am besten angepassten Abbildungsfunktion wiedergibt. Die derart 
ermittelten Korrekturwerte lassen sich in einer Korrekturtabelle speichern und fur eine 
Weiterverarbeitung nutzen. Es wird die am besten angepasste Abbildungsfunktion 
5 ermittelt. In vielen Fallen ist jedoch ein Restfehler der am besten angepassten 
Abbildungsfunktion gegenuber einer dem Abbildungsvorgang tatsachlich zugrunde 
liegenden Funktion unvermeidbar und/oder nur mit einem sehr hohen Aufwand 
vermeidbar. Mit Hilfe der aus der Korrekturtabelle abgeleiteten Korrekturfunktion ist 
eine Entzerrung mit hoher Qualitat auch dann moglich, wenn eine am besten 
1 0 angepasste Abbildungsfunktion einen Restfehler aufweist. 

In einer weiteren Ausfuhrungsform werden die Eckkoordinaten projizierter Bildelemente 
ermittelt. Eine Bildposition der Eckkoordinaten ergibt sich aus einer Projektion mit der 
ermittelten Abbildungsfunktion und der Korrekturwerte. Durch Interpolation der 
1 5 Korrekturwerte d er P asspunkte I asst s ich K orrekturfunktion a ufstellen. S omit konnen 
die Korrekturwerte an den Eckkoordinaten ermittelt werden. 

In einer weiteren Ausfuhrungsform wird ein Resampling derart durchgefuhrt, dass die 
Eckkoordinaten ein Vieleck, bevorzugt ein Viereck, aufspannen und die Grauwerte aller 
20 in dem Vieleck liegenden Bildelemente entsprechend ihrer GrolJe gewichtet gemittelt 
werden. 

In einer weiteren bevorzugten Ausfuhrungsform werden Objekte durch eine 
Klassifikation und/oder eine geometrische Strukturanalyse extrahiert. In der 

25 Klassifikation werden die Eigenschaften des Bildes analysiert und Objekte und/oder 
Flachen gleicher Klassenzugehorigkeit gebildet. Als Objekte und/oder Flachen gleicher 
Klassenzugehorigkeit werden Bildteile bezeichnet, deren Elemente (Pixel) jeweils einer 
Klasse angehoren. Die Klassifikation stutzt sich auf Bildeigenschaften. Diese 
Eigenschaften konnen sich auf einzelne Pixel beziehen. Eine Eigenschaft, die sich auf 

30 ein Pixel bezieht, ist beispielsweise ein Grauwert in einem Kanal, der durch das 
Reflexionsverhalten der dem Pixel entsprechenden Oberflache in einem festgelegten 
Wellenlangenbereich bestimmt wird. Weiter konnen Eigenschaften sich auch auf 
mehrere Pixel beziehen, bei denen die raumliche Anordnung einzelner Grauwerte eine 
Rolle spielt, beispielsweise bei der Textur. Eine weitere Eigenschaft ergibt sich aus 

35 dem Reflexionsverhalten der Oberflache in verschiedenen Wellenlangenbereichen 
(spektrale Eigenschaft). Ausgehend von bestimmten Eigenschaften lassen sich 
Objekte in Klassen einteilen, beispielsweise in Landnutzungsklassen wie Ackerbau, 
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Nadelwald, etc. Insbesondere werden damit Wolken und Gebiete mit stabilem 
Reflexionsverhalten, beispielsweise Wasser oder Nadelwalder, erkannt. Mittels einer 
Maske lassen sich die Bildobjekte, welche als Wolken erkannt wurden, kennzeichnen. 
Wolkengebiete werden fur die weitere Bearbeitung ausgeschlossen oder bei Bedarf 
5 einer gesonderten Bearbeitung unterworfen. Liegen multispektrale 
Fernerkundungsdaten vor, so werden Objekte vorzugsweise gemafi einem aus der 
Patentschrift DE 199 39 732 C2 bekannten Verfahren zur Klassifikation detektiert und 
identifiziert. 

10 Die Anhand ihrer Eigenschaften klassifizierten Objekte sind durch eine daran 
anschliefiende geometrische Strukturanalyse genauer klassifizierbar. Die geometrische 
Strukturanalyse umfasst beispielsweise eine Analyse der Objekte selbst, wobei deren 
Grofte, Form und/oder Homogenitat betrachtet werden. Vorzugsweise werden 
Informationen uber benachbarte Objekte einbezogen. So lassen sich beispielsweise 

15 Bildelemente in Ubergangsbereichen einem Objekt zuordnen. 

Hierzu werden im ersten Schritt die Randverlaufe der geeigneten klassifizierten 
Flachen ermittelt, wobei zusatzlich eine numerische Charakterisierung dieser Objekte 
durch ein geeignetes Strukturmali erfolgt. Dabei lassen sich die Objekte, welche an 
20 Wolken angrenzen, gesondert kennzeichnen. Diese Objekte werden dann mit 
entsprechenden Objekten des Referenzbildes korreliert, so dass als Ergebnis die 
Objekte vorab extrahiert und identifiziert sind. 

Aus den extrahierten Objekten werden die Positionen der jeweiligen abgeleiteten 
25 charakteristischen Objektpunkte, vorzugsweise Objektmittelpunkte, berechnet. 
Gegebenenfalls wird bei teilweise durch Wolken abgedeckten Objekten nach 
charakteristischen Randverlaufen gesucht. 

Die Bilder werden durch ein beliebiges Aufnahmegerat erzeugt, welches einkanalige 
30 oder mehrkanalige Bilder produziert. Ein geeignetes Aufnahmegerat ist beispielsweise 

eine Pushbroom-Camera. Das Aufnahmegerat ist an einem beliebigen Objekt 

'i 

angebracht, welches gegenuber einer aufzunehmenden Oberflache eine 
Relativbewegung d urchfuhrt. D abei i st e s a uch denkbar, d ass d ie O berflache relativ 
zum Aufnahmegerat bewegt wird. 

35 

Die erfindungsgemafie Vorrichtung umfasst einzelne Module. Die einzelnen Module 
der kpnnen entweder auf einem gemeinsamen Rechenmittel oder auf getrennten 
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Rechenmitteln integriert werden. Geeignete Rechenmittel sind beispielsweise ein PC, 
eine Workstation und/oder ein Rechnercluster. Die Verfahrensschritte konnen teilweise 
parallel bearbeitet werden. 

Die Erfindung wird nachfolgend anhand eines bevorzugten Ausfuhrungsbeispiels naher 
erlautert. Die Figuren zeigen: 

Fig. 1: den Ablauf bei einer automatischen Entzerrung eines Satellitenbildes, 
Fig. 2: eine schematische Darstellung der Generierung von Objektrandern 

durch Klassifikation und geometrische Strukturanalyse, 
Fig. 3: s chematisch einen Ausschnitt aus einem Bild sowie einem Referenzbild, 
Fig. 4: eine schematisch Illustration zur Ermittlung einer Abbildungsfunktion 

anhand geometrische Beziehungen, 
Fig. 5: e ine schematisch Illustration zur Ermittlung einer verbesserten 

Abbildungsfunktion anhand geometrische Beziehungen, 
Fig. 6: e ine schematische Darstellung von Korrekturwerten fur eine 

Abbildungsfunktion und 
Fig. 7: eine schematische Darstellung einer Entzerrung anhand einer 

Abbildungsfunktion. 

Fig. 1 zeigt schematisch den Ablauf bei einer automatischen Entzerrung eines 
einkanaligen oder mehrkanaligen Bildes auf ein Referenzbild. Das Verfahren ist 
modular durch die Module 1 bis 10 aufgebaut. Alle fur eine Entzerrung notwendigen 
Informationen sind aus Daten des Bildes und/oder eines Referenzbildes extrahierbar. 
Die Entzerrung auf das Referenzbild ist somit ohne Kenntnis der Abbildungsfunktion 
moglich. Liegen zusatzliche Informationen vor, so konnen diese im Verfahren 
verwendet werden. Fur eine Entzerrung ist es ausreichend, dass das Bild und das 
Referenzbild als einkanaliges Grauwertbild vorliegen. In den Modulen 1 bis 4 wird mit 
Hilfe von klassifizierten Objekten und deren Beziehungen zueinander eine 
Grobzuordnung zwischen Bild und Referenzbild erreicht. Auf der Grundlage der 
Grobzuordnung wird in Modul 5 ein grob angepasster virtueller Abbildungsvorgang 
ermittelt. In den Modulen 6 bis 8 erfolgt eine Feinanpassung des virtuellen 
Abbildungsvorganges. Als Ergebnis der Module 6 bis 8 liegt eine endgultige 
Abbildungsfunktion mit Korrekturtabelle vor. Die Zuordnung von Grauwertanteilen zu 
einer Flache (Pixel) findet in den Schritten 9 und 10 start. 
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Das Verfahren ist am Beispiel der Entzerrung eines Satellitenbildes S illustriert. Das 
Satellitenbild S wird beispielsweise durch eine Pushbroom-Camera aufgenommen. Die 
Pushbroom-camera besteht aus einer Sensorzeile. Durch eine Vorwartsbewegung des 
Satelliten Oder eines anderen bewegten Objekts beispielsweise eines Flugzeugs 
5 entsteht das Satellitenbild durch Aneinanderreihung der einzelnen Zeilen. 

Im Modul 1 wird zunachst das Satellitenbild S auf bestimmte Bildeigenschaften 
untersucht. Eigenschaften sind beispielsweise der Albedo-Wert, d.h. das Verhaltnis 
von reflektierter Strahlung zu Einstrahlung, Farb- oder Grauwerte und/oder texturelle 

10 Eigenschaften einer aufgenommenen Oberflache. Ausgehend von Eigenschaften 
lassen sich Bildobjekte in Klassen einteilen, beispielsweise in Landnutzungsklassen 
wie Ackerbau, Nadelwald, etc. Bildanteile werden somit im nachsten Schritt gemafi 
ihrer Eigenschaften klassifiziert und Flachen gleicher Klassenzugehorigkeit gebildet. 
Fur die Klassifikation wird ein geeigneter Klassifizierungsalgorithmus verwendet, durch 

15 welchen beispielsweise Gebiete mit stabilem Reflexionsverhalten, wie Wasserflachen, 
identifizierbar sind. 

In Modul 2 erfolgt eine Identifikation von Objekten, wobei Objektrander Randverlaufe 
der gebildeten Flachen sind und beispielsweise unter Verwendung eines 
20 Strukturanalysemodells eine numerische Charakterisierung eines Objekts aufgrund 
seines Strukturmafies erfolgt. Die Objektrandverlaufe ermbglichen eine 
Formbeschreibung der einzelnen Objekte. Die verschiedenen Objektformen sind eine 
Voraussetzung, urn gleiche Objekte in Referenzbild und Satellitenbild zu ermitteln. 

25 Fig. 2 illustriert schematised wie durch Klassifikation und geometrische 
Strukturanalyse ein Bild O ermittelt wird, dem einzelne Objektrander entnehmbar sind. 

In Modul 3 werden fur die aufgefundenen Objektrander Objektmittelpunkte berechnet, 
welche Position der Objekte wiedergeben. Das Satellitenbild S wird durch eine nicht 
30 Oder nur teilweise bekannte Abbildungsfunktion auf ein Referenzbild entzerrt. Die 
Objektmittelpunkte sowie weitere markante Punkte der Objektgeometrie sind 
Passpunkte fur eine Anpassung der Abbildungsfunktion. 

In dem Modul 4 werden markante extrahierte Objekte gemafi Bild O jeweils einem 
35 Objekt des Referenzbildes R zugeordnet und/oder extrahierte Objekte identifiziert. 
Dabei ist zu beachten, dass Objekte verzerrt dargestellt sind und daher eine 
Zuordnung aufgrund von Objektflachenformen nur bedingt moglich ist. Die Zuordnung 
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und/oder Identifikation erfolgt daher erfindungsgemafi zusatzlich durch Vektornetze, 
durch welche Objekte im Bild O sowie im Referenzbild R verknupft werden. Zum 
Erzeugen eines Vektornetzes sind mindestens drei markante Objekte notwendig. 
Durch eine Analyse der geometrischen Ahnlichkeit der Vektornetze werden Objekte 
5 raumlich zugeordnet und/oder identifiziert. 

Fig. 3 zeigt schematisch einen Ausschnitt aus einem Bild O sowie einem Referenzbild 
R. In dem Bild O sowie im Referenzbild R sind jeweils drei markante Objekten 01, 02, 
03 und Or, 02\ 03' aufgefunden, ihre Mittelpunkte als markante Objektpunkte 

10 identifiziert und als Passpunkte gewahlt. Die Passpunkte sind jeweils uber ein 
Vektornetz verbunden. Die Objekte 01, 02, 03 und das Vektornetz des Bildes O sind 
gegenuber dem Referenzbild verzerrt. Die Objekte 01, 02, 03 sind daher nicht direkt 
den Objekten 01\ 02', 03' zuordenbar und somit identifizierbar. Durch eine Analyse 
der geometrischen Ahnlichkeit der Vektornetze werden Objekte 01, 02, 03 raumlich 

1 5 den Objekten 01 ', 02\ 03' zugeordnet und dadurch identifiziert. 

In Modul 5 wird zunachst eine Anpassung der nicht oder nur teilweise bekannten 
Abbildungsfunktion an die in den Modulen 1 bis 4 ermittelten Passpunkte durchgefuhrt. 
Der Anfangsparametersatz wird dabei aus den geometrischen Beziehungen der 
20 Objekte zueinander berechnet. Die Anpassung erfolgt durch eine Anderung der 
Parameter der Abbildungsfunktion. 

Fig. 4 illustriert eine Bestimmung einer Abbildungsfunktion anhand geometrische 
Beziehungen. Bilder von der Erdoberflache 22 werden mittels eines Satelliten 20 
25 aufgenommen. Gesucht wird eine Abbildungsfunktion, durch welche die geometrischen 
Beziehungen von Objekten 01-04 der Erdoberflache zueinander moglichst exakt 
beibehalten werden 

In Modul 6 werden weitere Passpunkte fur eine Anpassung der Abbildungsfunktion 
30 generiert. Prinzipiell ist jeder Punkt, welcher sich in beiden Bildern eindeutig 
identifizieren lasst, als Passpunkt geeignet. Die potentiellen Passpunkte sind jedoch 
hinsichtlich ihrer Eignung unterschiedlich zu bewerten. Hierfur werden in Modul 6 
Qualitatskennzahlen von Passpunktstrukturen ermittelt. Passpunktstrukturen sind 
Flachen mit N*N Pixel, welche um einen potentiellen Passpunkt liegen. Die Grofie der 
35 Passpunktstrukturen ist ein Kompromiss zwischen dem idealen, dimensionslosen 
Pass-„Punkt" und einer wiedererkennbaren Bildstruktur. Geeignete 
Passpunktstrukturen sind beispielsweise 3x3 oder 5x5 Pixel. Die Qualitatskennzahlen 



11 



ergeben sich aus Veranderlichkeit, Richtungskontrast und Ahnlichkeit. Die 
Veranderlichkeit einer Passpunktstruktur wird aus ihrer Klassenzugehorigkeit 
entsprechend Modul 1 abgeleitet. Eine hohe Veranderlichkeit wirkt sich dabei negativ 
auf die Qualitatsbewertung des Passpunktes aus. Um in einem nachfolgenden 
5 Verfahrensschritt die Position des Passpunktes mit einer hohen Genauigkeit zu 
ermitteln, ist ein grofier Kontrast innerhalb der Passpunktstruktur von Vorteil. Einer 
Struktur mit grofien Richtungskontrasten wird daher ein hohes QualitatsmafJ 
zugeordnet. Die Ahnlichkeit zwischen zwei Strukturen ist fur einen Vergleich des 
projizierten Bildes mit dem Referenzbild sehr wichtig, da bei ungenugender Ahnlichkeit 
10 eine Positioniergenauigkeit stark absinkt. Das QualitatsmafJ ist daher proportional zur 
Ahnlichkeit. Jedem Passpunkt wird entsprechend seiner Qualitatskennzahl eine 
Wichtung zugeordnet. Passpunktstrukturen und dazugehorige Qualitatskennzahlen 
werden in einer Passpunkttabelle abgelegt. 

15 Mit der grofJeren Passpunktzahl und deren Wichtung wird in Modul 7 erneut eine 
Anpassung der Parameter der Abbildungsfunktion durchgefuhrt, wodurch die 
Genauigkeit der Abbildungsfunktion verbessert wird. 

Fig. 5 illustriert eine Wirkung der Anpassung von Parametern fur eine verbesserte 
20 Bestimmung der Abbildungsfunktion. Die Abbildungsfunktion wird durch eine Vielzahl 
an Gegebenheiten zum Zeitpunkt einer Aufnahme bestimmt. So hat beispielsweise 
eine Ausrichtung des Satelliten 20 eine Auswirkung auf eine Abbildung der 
Erdoberflache 22. Die Ausrichtung des Satelliten 20 sowie andere Gegebenheiten sind 
jedoch im Regelfall nur in Grenzen bekannt. Die Abbildungsfunktion wird durch eine 
25 Variation mindestens eines Parameters so angepasst, dass die Positionsfehler im Bild 
minimiert werden. Die Anpassung erfolgt, ohne dass die Ausrichtung des Satelliten 20 
und/oder andere Gegebenheiten zum Zeitpunkt der Aufnahme bekannt sind. 

Positionsrestfehlern, welche durch die verbesserte Abbildungsfunktion nicht vollstandig 
30 vermieden werden, wird mit einer Ausgleichsrechnung in Modul 8 begegnet. Dabei wird 
angenommen, dass die Positionsrestfehler uber kurze Distanzen keine starke 
Anderung aufweisen. Die Ausgleichsrechnung wird daher unter Einbeziehung der 
naheren Umgebung durchgefuhrt. Als Ergebnis der Ausgleichsrechnung wird in Modul 
8 eine Fehlerkorrekturtabelle erstellt. Durch Interpolation der Fehlerkorrekturtabelle ist 
35 eine Fehlerkorrekturfunktion ermittelbar. 
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Fig. 6 verdeutlicht die Korrekturtabelle. Die Korrekturwerte sind durch Pfeile 
dargestellt. Einzelnen Objekten 06 - 08 sind unterschiedliche Korrekturwerte 
zugeordnet, was durch unterschiedliche Pfeillange illustriert wird. 

In Modul 9 werden Eckkoordinaten der projizierten Bildelemente berechnet. Die 
Berechnung erfolgt unter Verwendung der ermittelten Abbildungsfunktion und der 
Fehlerkorrekturfunktion. Zum Schluss findet in Modul 10 ein Resampling mit den vorher 
berechneten projizierten Positionen statt. 

Fig. 7 illustriert eine Entzerrung des Bildes anhand der ermittelten Abbildungsfunktion 
mit Fehlerkorrektur, wobei eine Bildstruktur B1 auf eine Referenzstruktur R1 projiziert 
wird. 



PATENTANSPRUCHE 



Verfahren zur automatischen Entzerrung von Bildern, wobei mindestens ein 
Bild durch eine Abbildungsfunktion auf ein Referenzbild (R) entzerrt wird und 
Parameter der Abbildungsfunktion mindestens teilweise unbekannt sind, 
umfassend mindestens 

eine Extraktion von mindestens drei Objekten (01 - 03) aus dem Bild 

(O), 

eine Ermittlung von mindestens drei Passpunkten in dem Bild, wobei 
charakteristische Objektpunkte der extrahierten Objekte als Passpunkte 
ermittelt werden, 

- eine Zuordnung der Objekte (01 - 03) zu Objekten (Or - 03') im 
Referenzbild, wobei eine Zuordnung von Objekten in beiden Bildern 
aufgrund einer Ahnlichkeit von Objekten und/oder eines Vektornetzes 
erfolgt und das Vektornetz durch eine Verbindung charakteristischer 
Objektpunkte gebildet wird, und 

eine Auswahl einer geeigneten Abbildungsfunktion und/oder eine 
Anpassung der Parameter der Abbildungsfunktion, wobei die 
Abbildungsfunktion durch eine Anderung der Parameter derart verandert 
wird, dass ein Fehlersummenminimum in Positionsdifferenzen zwischen 
projizierten Passpunkten und entsprechenden Punkten im Referenzbild 
erreicht wird. 

Verfahren nach Anspruch 1, umfassend eine Generierung von gewichteten 
Passpunkten, wobei urn einen Passpunkt des Bildes und/oder des 
Referenzbildes eine Passpunktstruktur umfassend eine begrenzte Anzahl an 
Pixeln gebildet wird, die Passpunktstruktur durch die Abbildungsfunktion auf 
das jeweils andere Bild als Bildstruktur projiziert wird, um zu erkennen, ob auch 
dort eine entsprechende Bildstruktur ausreichender Qualitat vorhanden ist, 
wobei eine Qualitat der Passpunktstruktur mindestens durch Veranderlichkeit, 
Richtungskontrast und/oder Ahnlichkeit beschrieben ist und aufgrund der 
Qualitat eine Wichtung der Passpunkte ermittelt wird. 

Verfahren nach Anspruch 2, umfassend eine Anpassung der Position des 
Passpunktes im Bild und/oder im Referenzbild, wobei fur mindestens einen 
Kanal die Form einer Grauwertverteilung der Passpunktstruktur im Referenzbild 
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und die Form einer Grauwertverteilung der Bildstruktur im Bild aufeinander 
abgestimmt werden, im Bild und/oder im Referenzbild mindestens eine 
Differenz zwischen Grauwerten zweier benachbarter Pixel der 
Passpunktstruktur und mindestens eine Differenz zwischen Grauwerten der 
entsprechenden Pixel der Bildstruktur gebildet wird und aus einem Unterschied 
der Differenzen ein Fehler abgeleitet wird, das unscharfe Strukturteil in das 
scharfe Strukturteil abgebildet wird, die Passpunktstruktur im Bild und/oder im 
Referenzbild in vertikaler und/oder horizontaler Richtung verschoben wird und 
fur die neue Position der Fehler ermittelt wird. 

Verfahren nach Anspruch 2 oder 3, umfassend eine Anpassung der Parameter 
der Abbildungsfunktion und/oder eine Auswahl einer geeigneten 
Abbildungsfunktion, wobei die Abbildungsfunktion durch eine Anderung der 
Parameter derart verandert wird, dass ein Fehlersummenminimum in 
Positionsdifferenzen zwischen projizierten Passpunkten und entsprechenden 
Passpunkten im Referenzbild unter Berucksichtigung der Wichtung der 
Passpunkte erreicht wird. 

Verfahren nach einem der Anspruche 1 bis 4, umfassend eine 
Ausgleichsrechnung mittels einer Korrekturfunktion, wobei fur mindestens zwei 
Passpunkte mindestens ein vertikaler und/oder ein horizontaler Korrekturwert 
bestimmt werden, welche die Positionsdifferenz zwischen dem projizierten 
Passpunkt und dem entsprechenden Passpunkt im Referenzbild angeben, und 
die Korrekturfunktion in Abhangigkeit der Korrekturwerte bestimmt wird. 

Verfahren nach Anspruch 5, umfassend eine Projektion von Eckkoordinaten 
eines Bildelementes auf Bildpositionen, wobei die Bildpositionen der 
Eckkoordinaten aufgrund der Abbildungsfunktion und der Korrekturfunktion 
ermittelt werden. 

Verfahren nach Anspruch 6, umfassend ein Resampling, wobei die 
Eckkoordinaten ein Vieleck, bevorzugt ein Viereck, aufspannen und die 
Grauwerte entsprechend der Flachenanteile aller in dem Vieleck liegenden 
Bildelemente in den resultierenden Grauwert eingehen. 



. Verfahren nach einem der Anspruche 1 bis 7, dadurch gekennzeichnet, dass 
die Extraktion eine Klassifikation und/oder eine geometrische Strukturanalyse 
umfasst, wobei 

in der Klassifikation die Eigenschaften des Bildes analysiert werden und 
Objekte und/oder Flachen gleicher Klassenzugehorigkeit gebildet 
werden und 

in der geometrischen Strukturanalyse ein Randverlauf eines Objekts aus 
einem Randverlauf einer Flache ermittelt wird und/oder eine numerische 
Charakterisierung der Objekte durch ein StrukturmafJ erfolgt. 

. Vorrichtung zur automatischen Entzerrung von Bildern, wobei mindestens ein 
Bild durch eine Abbildungsfunktion auf ein Referenzbild (R) entzerrbar ist und 
Parameter der Abbildungsfunktion mindestens teilweise unbekannt sind, 
umfassend mindestens 

- ein Modul (1, 2) fur eine Extraktion von mindestens drei Objekten (01 - 
03) aus dem Bild (O), 

- ein Modul (3) fur eine Ermittlung von mindestens drei Passpunkten in 
dem Bild, wobei charakteristische Objektpunkte der extrahierten Objekte 
als Passpunkte ermittelbar sind, 

- ein Modul (4) fur eine Zuordnung der Objekte (01 - 03) zu Objekten 
(OV - 03') im Referenzbild, wobei eine Zuordnung von Objekten in 
beiden Bildern aufgrund einer Ahnlichkeit von Objekten und/oder eines 
Vektornetzes erfolgt und das Vektornetz durch eine Verbindung 
charakteristischer Objektpunkte bildbar ist, und 

- eine Auswahl einer geeigneten Abbildungsfunktion und/oder eine 
Anpassung der Parameter der Abbildungsfunktion, wobei die 
Abbildungsfunktion durch eine Anderung der Parameter derart verandert 
wird, dass ein Fehlersummenminimum in Positionsdifferenzen zwischen 
projizierten Passpunkten und entsprechenden Punkten im Referenzbild 
erreicht wird. 

0. Vorrichtung nach Anspruch 9, umfassend ein Modul (6) fur eine Generierung 
gewichteter Passpunkte, durch welches urn einen Passpunkt des Bildes 
und/oder des Referenzbildes eine Passpunktstruktur umfassend eine 
begrenzte Anzahl an Pixel bildbar ist, die Passpunktstruktur durch die 
Abbildungsfunktion auf das jeweils andere Bild als Bildstruktur abbildbar ist, 
wobei eine Qualitat der Passpunktstruktur mindestens durch Veranderlichkeit, 
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Richtungskontrast und/oder Ahnlichkeit beschreibbar ist und aufgrund der 
Qualitat eine Wichtung der Passpunkte ermittelbar ist. 

1 1 . Vorrichtung nach Anspruch 10, umfassend ein Modul, durch welches eine 
Position des Passpunktes im Bild und/oder im Referenzbild anpassbar ist, 
wobei auf mindestens einem Kanal die Form einer Grauwertverteilung der 
Passpunktstruktur und der Bildstruktur aufeinander abstimmbar sind, 
mindestens eine Differenz zwischen Grauwerten zweier benachbarter Pixel der 
Passpunktstruktur und mindestens eine Differenz zwischen Grauwerten der 
entsprechenden Pixel der Bildstruktur bildbar ist und aus einem Unterschied der 
Differenzen ein Fehler ableitbar ist, wobei das unscharfe Strukturteil in das 
scharfe Strukturteil abgebildet wird, die Passpunktstruktur im Bild und/oder im 
Referenzbild in vertikaler und/oder horizontaler Richtung verschiebbar ist und 
fur die neue Position der Fehler ermittelbar ist. 

12. Vorrichtung nach Anspruch 10 oder 11, umfassend ein Modul (7), durch 
welches die Parameter der Abbildungsfunktion anpassbar sind, wobei. die 
Abbildungsfunktion durch eine Anderung der Parameter derart veranderbar ist, 
dass ein Fehlersummenminimum in den Positionsdifferenzen zwischen 
Passpunkten und zugehorigen projizierten Bildpunkten unter Berucksichtigung 
der Wichtung der Passpunkte minimiert wird. 

13. Vorrichtung nach einem der Anspruche 9 bis 12, umfassend ein Modul (8), 
durch welches eine Ausgleichsrechnung durchfuhrbar ist, wobei fur jeden 
Passpunkt mindestens ein vertikaler und/oder ein horizontaler Korrekturwert 
bestimmbar ist, welcher die Abweichung des Wertes der Abbildungsfunktion 
vom Wert der Ausgleichsfunktion am Ort des Passpunktes wiedergibt. 

14. Vorrichtung nach einem der Anspruche 9 bis 13, umfassend ein Modul (9), 
durch welches Eckkoordinaten eines Bildelementes auf Bildpositionen abbildbar 
sind, wobei die Bildpositionen der Eckkoordinaten aufgrund der 
Abbildungsfunktion und der Korrekturfunktion ermittelbar sind. 

1 5. Vorrichtung n ach e inem d er A nspruche 9 b is 1 4, u mfassend e in M odul ( 1 0), 
durch welches ein Resampling durchfuhrbar ist, wobei die Eckkoordinaten ein 
Vieleck, b evorzugt e in Viereck, a ufspannen u nd d ie G rauwerte e ntsprechend 
der Flachenanteile aller in dem Vieleck liegenden Bildelemente den 
resultierenden Grauwert bestimmen. 
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16. Vorrichtung nach einem der Anspruche 9 bis 15, dadurch gekennzeichnet, dass 
die Extraktion ein Modul (1) fur eine Klassifikation und/oder ein Modul (2) fur 
eine geometrische Strukturanalyse umfasst, wobei 

in der Klassifikation die Eigenschaften des Bildes analysierbar sind und 
Objekte und/oder Flachen gleicher Klassenzugehorigkeit bildbar und 
in der geometrischen Strukturanalyse ein Randverlauf eines Objekts aus 
einem Randverlauf einer Flache ermittelbar ist und/oder eine 
numerische Charakterisierung eines Objekts durch ein Strukturmafi 
erfolgt. 
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ZUSAMMENFASSUNG 



Verfahren und Vorrichtung zur automatischen 
Entzerrung von Bildern 



Die Erfindung betrifft ein Verfahren und eine Vorrichtung automatischen Entzerrung 
von Bildern, wobei mindestens ein Bild durch eine Abbildungsfunktion auf ein 
Referenzbild entzerrt wird und Parameter der Abbildungsfunktion mindestens teilweise 
unbekannt sind, umfassend mindestens eine Extraktion von mindestens drei Objekten 
(01 - 03) aus dem Bild (O), eine Ermittlung von mindestens drei Passpunkten in dem 
Bild, wobei charakteristische Objektpunkte der extrahierten Objekte als Passpunkte 
ermittelt werden, eine Zuordnung der Objekte (01 - 03) zu Objekten (OV - 03') im 
Referenzbild, wobei eine Zuordnung von Objekten in beiden Bildern aufgrund einer 
Ahnlichkeit von Objekten und/oder eines Vektornetzes erfolgt und das Vektornetz 
durch eine Verbindung charakteristischer Objektpunkte gebildet wird, und eine 
Auswahl einer geeigneten Abbildungsfunktion und/oder eine Anpassung der Parameter 
der Abbildungsfunktion, wobei die Abbildungsfunktion durch eine Anderung der 
Parameter derart verandert wird, dass ein Fehlersummenminimum in 
Positionsdifferenzen zwischen projizierten Passpunkten und entsprechenden Punkten 
im Referenzbild erreicht wird. 



(Fig. 1) 



